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Abstract 

The mass of a sphere of simmetry in the Lemaitre universe is discussed using 
the Hawking-Hay ward quasi-local energy and clarifying existing ambiguities. A 
covariantly conserved current introduced by Cahill and McVittie is shown to be a 
multiple of the Kodama energy current. 
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1 Introduction 


The Lemaitre model [lj is a spherically symmetric inhomogeneous universe which solves 
the Einstein equations and generalizes the better known Lemaitre-Tolman-Bondi (LTB) 
|2[ [3] geometry to the case of non-zero pressure. Inhomogeneous universes (see Ref. [0J 
for a comprehensive review) have been the subject of much work in recent years because 
of the attempts to explain the current acceleration of the universe without an ad hoc 
dark energy or cosmological constant and without abandoning general relativity (see [5] 
for a review). This work on LTB models has revived interest also in Lemaitre models 
and inhomogeneous cosmologies in general. The issue has been raised recently of what 
is the physical mass contained in a sphere in the Lemaitre universe [6]. Over the years, 
there have been various proposals (reviewed in [6]). More recently, LTB models have 
been studied in various contexts [7] and the study has moved to a higher degree of 
sophistication with the introduction of perturbations of LTB models [8] ; similar to the 
case of Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetimes, perturbations allow 
research on more sophisticated physics and open the door to the use of tools such as the 
temperature fluctuations in the cosmic microwave background [9]. As a starting point 
for this and future developments, it would be reassuring to know that we understand the 
physics of the unperturbed Lemaitre universe, including its energy. Here we show that 
approaching this issue using the Hawking-Hayward quasi-local energy US HU and the 
Kodarna vector |T2] dissipates some ambiguities in the recent literature. In the presence 
of spherical symmetry (and, therefore, also in the Lemaitre universe), the Hawking- 
Hayward quasi-local energy reduces HUE! to the Misner-Sharp-Hernandez mass [T3] . 
Spherical symmetry allows one to introduce also the Kodarna vector [T2j, which generates 
a covariantly conserved current and has the Misner-Sharp-Hernandez mass as its Noether 
charge nu- 

in the next section we discuss the Lemaitre universe and the Misner-Sharp-Hernandez 
mass contained in sphere of symmetry^ showing how this clarifies doubts that still linger 
on in the literature. The following section discusses the Kodarna vector and shows that 
a conserved current discussed in Refs. HUE] is a multiple of the Kodarna current built 
out of the Einstein tensor and the Kodarna vector. 

We use units in which the speed of light and Newton’s constant are unity and we 
follow the notation of Wald’s textbook [15] . 

: By “sphere of symmetry” we mean a 2-dimensional surface which is an orbit of the isometry of the 
spacetime manifold describing spherical symmetry (of course, such orbits exist through any point of a 
spherically symmetric spacetime). 
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2 Lemaitre geometry and the mass of a sphere of 
symmetry 

The spherically symmetric Lemaitre line element in coordinates {aU} = {t, r, 9, ip} co¬ 
moving with the fluid source is [I] 

ds 2 = —e 2a dt 2 + e x dr 2 + R 2 dQ 2 2 ^, (2.1) 

where a = a(t,r), A = A(t,r), and R(t,r) is the areal radius, while dQ 2 ^ = dO 2 + 
sin 2 9 dip 2 is the line element on the unit 2-sphere. The simplest stress-energy tensor 
sourcing the Lemaitre spacetime is that of a perfect fluid 

T ab = (P + p) u a u b + Pg ab (2.2) 

where p(t, r) and P(t, r ) are the energy density and pressure of the fluid as perceived by 
a comoving observer with 4-velocity u c with components u M = (e _<7 , 0, 0, 0) in comoving 
coordinates. In general, however, the matter source of the Lemaitre metric (12.11) is not 
restricted to be a perfect fluid. 

Following [6], we assume isotropic pressure (which is a restriction on the full Lemaitre 
spacetime) and we allow for a cosmological constant A. The Lemaitre universe is a 
solution of the Einstein equations 


(Tab T A g ab ^^Tah , 


(2.3) 


where G ab is the Einstein tensor and g ab is the spacetime metric. 

In order to discuss the Misner-Sharp-Hernandez mass Mmsh and the Kodama vector 
k a , it is convenient to express the line element (12.ip using the areal radius R = x 1 
instead of r = x 1 because then Mmsh and k a assume straightforward expressions. Our 
conclusions, however, are fully coordinate-independent and it will be easy to revert to 
the coordinates {M 1 } = {£, r, 9,<p}. Our goal is to recast the line element (12. Tj) in the 
form 

ds 2 = —A(T, R)dT 2 + B(T, R)dR 2 + R 2 dn 2 {2) . (2.4) 

Using the differential relation dR = Rdt + R'dr, where an overdot and a prime denote 
differentiation with respect to t and r, respectively, the line element (12.1|) becomes 


ds 2 



dt 2 + —dR 2 
R!- 


2 Re x 
R’ 2 


dtdR + R 2 dfl 2 2 ^. 


(2.5) 
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The cross-term in dtclR is eliminated by introducing a new time coordinate x 0 ' = T 
defined by 

dT = 4 (,dt + pdR) , (2.6) 

r 

where P(t,r ) is a function to be determined and F(t,r) is an integrating factor which 
must be introduced to guarantee that dT is an exact differential. It satisfies the equation 


d (\ 


d (P 


dr \FJ dt \F / 

Upon substitution of dt = FdT — PdR, the line element (12.5ft assumes the form 


(2.7) 


ds‘ = - I e 2 ' - e A 


e A 

R 2 


R 2 
R ~ 2 


F 2 dT 2 + 2 F 


P \e 2a - e A 


R 2 




R’ 2 / R’ 2 


dTdR 


Jla ..A 


R 


R 


e A — p 2 + 2P—e 


Setting 

one obtains 

ds 2 = 


R' 2 


P(t,r) = 


R' 2 


dR 2 + R 2 dfl 2 


( 2 ) ■ 


Re* 


R' 2 ( e 2f7 - e A 


- e A ^ ) F 2 dT 2 + 


e A e 2cr 


R' 2 ( e 2fJ - e A §2 


dR 2 + R 2 dn 2 {2) , 


which is of the form (j2.4|) with 

A(T,R) = 

B(T,R ) = 


e 2a - e A 


R 2 

If 2 


02(j 


(.: 2 " - (.>■ y. 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 

( 2 . 12 ) 


In his original paper, Lemaitre identified (without a real justification) the mass contained 
in a sphere of radius r with the quantity 


R 


M = — I 1 + R 2 e~ 2a - i?'V A - 


Ai? 2 


(2.13) 
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The derivatives of this “mass” M are related to the cosmic fluid density and pressure 
by PEI 

M' = 4ttR 2 R' P , (2.14) 

M = -4? tR 2 RP. (2.15) 

In 1970, Cahill and McVittie [T4] studied the Lemaitre model with A = 0 and identified 
the right hand side of eq. (j2.13(1 (with A = 0) with the physical mass. The rationale was 
that, if a Lemaitre model is joined to an exterior Schwarschild geometry in a Swiss-cheese 
model, then M must equal the exterior massj§ and that the Bianchi identities generalize 
eqs. (12. 14j) and 02.151) . Cahill and McVittie were aware that their mass proposal M 
coincided with the Misner-Sharp-Hernandez mass, then just recently introduced [ITj . 

The authors of Ref. [6] note that the interpretation of M is rather clear when there is 
no pressure and A is absent, but it is not so straighforward otherwise. They proceed to 
note that eq. 02.141) . which does not contain P, holds in all Lemaitre models. This is a 
good observation but, in principle, there could be other quantities or equations which do 
not depend explicitly on P and could be used to define effective masses. As a guideline 
to find out the physical mass of sphere, the authors of |B] proceed to analyze the geodesic 
deviation equation in order to establish which mass is “seen” by test particles, and the 
discussion necessarily becomes rather involved. They reach the rather disheartening 
conclusion that “we cannot separate the mass, the cosmological constant, the density 
and the pressure from each other, and so we cannot create a unique definition of mass 
based on geometric invariants of the metric in the general case” [6]. Indeed, we can. 
Our knowledge of energy and mass in general relativity has progressed greatly since 
the times of Lemaitre and Cahill and McVittie, with the introduction of the various 
quasi-local energies (see [T6] for a review), which culminated in the Hawking-Hayward 
quasi-local energy mm- In spherical symmetry, the Hawking-Hayward quasi-local 
energy reduces hd ini to the Misner-Sharp-Hernandez mass Mmsh [E>j, which is defined 
in a coordinate-invariant way by Ha hh m m mm 

1 - 2M ^ SH = V C PV C P, (2.16) 

R 

where R is the areal radius (which, being related to the area A of 2-spheres of symmetry 
by R = is a geometrically defined quantity). 

2 In retrospect, this is a good argument because there is little arguing on the physical mass of the 
Schwarzschild spacetime, and the choice proved to give the correct answer (see below). 
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In coordinates { x 11 '} = {T, R,6,<p}, the squared gradient in the right hand side of 
eq. (I2.16p is simply g RR = B~ 1 and gives the Misner-Sharp-Hernandez mass 


Mmsh = 9 RR ) 

= ^ ^1 + R 2 e~ 2a — R' 2 e_A ) 

A R 3 

= M -\ -— • (2.17) 

o 


So, for A = 0, the Cahill-McVittie prescription coincides with the Hawking-Hayward/Misner- 
Sharp-Hernandez one. This result can, of course, be obtained also in coordinates 
{t, r, 9, <p} by computing V C KV C R with the metric (12.ip and using the well known rela¬ 
tion 


R 


±e° 


2 M 


R 


+ R' 2 e~ A - 1 + 


A R 2 


(2.18) 


which follows from the definition (12.131) [6j. 

Note that the Misner-Sharp-Hernandez mass does not depend explicitly on the pres¬ 
sure P (although P affects the cosmic expansion and determines the metric coefficients 
which, in turn, determine Mmsh), while M depends on P but not on p (eq. (12.15P ). 
This fact is well known in FLRW space mm, to which the Lemaitre model reduces 
if a = 1 and A = A (t). Therefore, there is no issue of disentangling the contribution 
of P from those of p and A. The contribution AR 3 /6 to Mmsh is easily interpreted as 
the mass corresponding to the cosmological constant energy density p\ = -^ contained 
in a sphere of areal radius R and volume AnR 3 /3. The decomposition of Mmsh hAo a 
“local” and a “cosmological” part is covariant: this point has been discussed in detail in 
[19) for the McVittie metric ra. which exhibits some of the properties of the Lemaitre 
model (although it belongs to a different family of solutions of the Einstein equations), 
and we will not repeat the discussion here. 

Although the Hawking-Hayward mass is not mentioned in [6j , the authors some¬ 
how end up reasoning along the same lines in their quest for the physical mass, when 
they stress the role of the apparent horizon in relating cosmic mass and diameter dis¬ 
tance. In fact, the apparent horizons of any spherically symmetric metric are defined by 
V C RV C R = 0 (e.g., [ 21] [22]). which relates the apparent horizon radii with the Misner- 
Sharp-Hernandez mass contained through the relation Rah = 2Mmsh (which mimics the 
expression of the Schwarzschild radius in the Schwarzschild spacetime) [21, 22]. In the 
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Lemaitre model, the recipe to locate the apparent horizons translates into g RR = 0, or 

i? ,2 e 2cr — R 2 e x = 0 . (2.19) 

In general, multiple solutions to this equation (which must be solved numerically) are 
possible, describing both black hole and cosmological apparent horizons which evolve in 
time (examples are solved in [[251123] ). 

Using the Misner-Sharp-Hernandez mass (12.17p . the line element (12.10)1 can now be 
written as 

* 2 = - e " +A (0 2 O'- 2J jr) dT2 + 0Sz5T + «*<«&> ■ (2.20) 

The spatial part of this line element resembles the spatial part of the Schwarzschild 
metric, but only superficially because Mmsh is not a constant but depends on the areal 
radius R. 


3 The Kodama energy current 

In a generic Lemaitre model there is no timelike Killing vector but, as in any spherically 
symmetric spacetime, one can introduce the closest thing to it, the Kodama vector na 

k‘ = e ab V b R, (3.1) 

where e ab is the volume element on the 2-space orthogonal to the 2-spheres of symmetry. 
If the spacetime metric is decomposed according to 

ds 2 = h ab dx a dx b + R 2 dQ 2 ^ (a, 6 = 0,1), (3.2) 

then e ab is the volume form associated with the 2-metric h ab . The Kodama vector plays 
the role of a Killing vector where there is none: it is timelike in a spacetime region, null 
on an apparent horizon, and becomes spacelike on the other side of this horizon ra- 
What makes the Kodama vector remarkable is the fact that the Kodama current 

J a = G ab k b (3.3) 

associated with it is covariantly conserved [12], V C J C = 0, such an unexpected property 
to be called the “Kodama miracle” |22j . What is more, the Misner-Sharp-Hernandez 
mass almost universally identified with the physical mass-energy in spherical symmetry 
in general relativity, turns out to be the Noether charge associated with the Kodama 
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current HU- There are strong claims in the literature that, in the absence of a preferred 
time derived from a timelike Killing vector, the Kodama vector introduces a preferred 
time and surface gravity on apparent horizons, which determine a Hawking temperature 
and make thermodynamics meaningful for time-evolving apparent horizons (see Ref. {25 , 
[261 27] for reviews). Given the geometry (12.4(1 . the Kodama vector is pT2] 


P 


-1 

\[AB 



(3.4) 


and, in the Lemaitre geometry (12.101) . it has components 

i^'le-^e-V 2 


P = 


F 


P 1 


(3.5) 


Cahill and McVittie found a conserved 4-current (reported also in [6]) which, in 
comoving coordinates {P} = has components 


_ sind 
H CM ) 4vr^ 


0,0) = 


-trp-A/2 


e e 


4vr R 2 


M\ —M, 0, 0 


(3.6) 


Eqs. (12.14ft and (12.15jl then yield 


J (CM) = e~ ff e ~ x/2 RP, 0, o) . (3.7) 

It is natural to ask whether this current is the same as the Kodama current. To find out, 
one computes the Kodama vector (13.11) . which is found to have components in comoving 
coordinates 

P = <rv v R = e^R + e^R'. (3.8) 

Therefore, it is 

k° = e 01 R r = g 00 g u \f\h\ R' = —e _f7 e _A/2 i?', (3.9) 

k 1 = e 10 R = -g m g ll -\J\h\R = e _ff e _A/2 R, (3.10) 

where h = —e 2cr e A is the determinant of the 2-metric h a b in the submanifold orthogonal 
to the 2-spheres of symmetry. We have 


p = e - CT e - A/2 R, 0 , o) 


(3.11) 
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and, lowering the indices, 


k n = (e a e- x/2 R',e~ a e x/2 R, 0,o) . (3.12) 

Using the Einstein equations (I2.3jl and the fluid four-velocity u M = — e CT S tl0l the non¬ 
vanishing components of the Kodarna current J a = G ab kb are found to be 

J° = G 00 k 0 = ( g 00 ) 2 G 00 k 0 = e~ a e~ x/2 R' (8vr p + A) , (3.13) 

J 1 = G u h = (g 11 ) 2 Gi 1 k 1 = e~ a e~ x/2 R (8vr P - A) , (3.14) 

By comparing the expression 

= e - ff e _A/2 (# (8vr p + A ),R (8?rP - A), 0, o) (3.15) 

with eq. (13.71) . it is clear that the Cahill-McVittie conserved current is just 

J(CU, = (Stt)- 1 ./' 1 (3.16) 

when A = 0 (and it is obtained by contracting T ab , instead of G ab , with the Kodarna 
vector). This fact was to be suspected, since asking for two separate “miracles” would 
be asking too much. 

4 Conclusions 

There are several quasi-local energy constructs in general relativity (see the review [16]) 
and there is no mathematical “proof” selecting the “correct” one. However, there is 
now general consensus that the Hawking-Hayward quasi-local construction ma na is 
preferred in the sense that it encapsulates better than its competitors the physical 
properties required by the mass-energy of a system. Among the advantages of the 
Hawking-Hayward mass notion are the facts that it is well defined for non-asymptotically 
flat and non-stationary spacetimes [T6|. In spherical symmetry, the Hawking-Hayward 
quasi-local energy reduces na to the better known Misner-Sharp-Hernandez mass [15] . 
It appears that Cahill and McVittie identified the correct mass notion in the A = 0 
Lemaitre space, the Hawking-Hayward/Misner-Sharp-Hernandez one. It is well known 
since the early days of the Hawking-Hayward mass mm that this object is also the 
Noether charge associated with the Kodarna current, and we won’t repeat the derivation 
of this result here. The Cahill-McVittie covariantly conserved current is just a multiple 
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of the Kodama energy current. Lemaitre, Cahill, and McVittie made a clever choice 
and they were correct after all, although perhaps they could not be sure of the reason 
why, because they came decades before Hawking and Hayward or before the present-day 
understanding of the quasi-local energy. 


Acknowledgments 

This work is supported by Bishop’s University and by the Natural Sciences and Engi¬ 
neering Research Council of Canada ( NSERC ). 


References 

[1] Lemaitre, G., Ann. Soc. Sci. Bruxelles A 53, 51 (1933) [reprinted in Gen. Rel. Grav. 
29, 641 (1997)] 

[2] Tolrnan, R.C., Proc. Nat. Acad. Sci U.S.A. 20, 169 (1934) 

[3] Bondi, H., Mon. Not. R. Astron. Soc. 107, 410 (1947) 

[4] Krasinski, A., Inhomogeneous Cosmological Models (Cambridge University Press, 
Cambridge, 1997) 

[5] Bolejko, K., Celerier, M.-N., and Krasinski, A., Class. Quantum Grav. 28, 164002 

( 2011 ) 

[6] Alfedeel, A.A.H. and Hcllaby, C., Gen. Rel. Grav. 42, 1935 (2010) 

[7] Pavlidou, V., Tetradis, N., and Tomaras, T.N., JCAP 1405, 017 (2014); Bolejko, 
K., Hcllaby, C., and Alfedeel, A.H.A., JCAP 1109, 011 (2011); Nogueira, F.A.M.G., 
arXiv: 1312.5005; Iribarrem, A. et al., arXiv: 1308.2199; Marra, V. and Paakkonen, 
M. JCAP 1201, 025 (2012); Clarkson, C. and Regis, M., JCAP 1102, 013 (2011); 
Clarkson, C. and Maartens, R., Class. Quantum Grav. 27, 124008 (2010); Moradi, 
R., Firouzjaee, J.T., and Mansouri, R., arXiv:1301.1480 

[8] Leithes, A. and Malik, K.A., arXiv: 1403.7661 

[9] Durrer, R., The Cosmic Microwave Background (Cambridge University Press, Cam¬ 
bridge, 2008) 

[10] Hawking, S.W., J. Math. Phys. 9, 598 (1968) 


9 


[11] Hayward, S.A., Phys. Rev. D 49 , 831 (1994) 

[12] Kodama, H., Prog. Theor. Phys. 63, 1217 (1980) 

[13] Misner, C.W. and Sharp, D.H., Phys. Rev. 136 , B571 (1964); Hernandez, W.C. and 
Misner, C.W., Astrophys. J. 143 , 452 (1966) 

[14] Cahill, M.E. and McVittie, G.C., J. Math. Phys. 11, 1382 (1970) 

[15] Wald, R.M., General Relativity (University of Chicago Press, Chicago, 1984) 

[16] Szabados, L. Living Rev. Relat. 7, 4 (2004) 

[17] Hayward, S.A., Phys. Rev. D 53, 1938 (1996) 

[18] Faraoni, V. Phys. Rev. D 84, 024003 (2011) 

[19] Carrera, M. and Giulini, D., Rev. Mod. Phys. 82, 169 (2010) 

[20] McVittie, G.C., Mon. Not. R. Astron. Soc. 93, 325 (1933) 

[21] Nielsen A.B. and Visser, Class. Quantum Grav. 23, 4637 (2006) 

[22] Abreu, G. and Visser, M., Phys. Rev. D 82 , 044027 (2010) 

[23] Booth, I., Brits, L, Gonzalez, J.A., and van den Broeck, V., Class. Quantum Grav. 
23 , 413 (2006) 

[24] Gao, C., Chen, X., Shen, Y.-G., and Faraoni, V., Phys. Rev. D 84 , 104047 (2011) 

[25] Di Criscienzo, R., Hayward, S.A., Nadalini, M., Vanzo, L., and Zerbini, S., Class. 
Quantum Grav. 27 , 015006 (2010) 

[26] Faraoni, V., Galaxies 1, 114 (2013) [arXiv: 1309.4915] 

[27] Faraoni, V. Cosmological and Black Hole Apparent Horizons (Springer, New York, 
forthcoming 2015) 


10 


